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ABSTRACT 

Different determinations currently suggest scattered values for the power spectrum 
normalization on the scale of galaxy clusters, as- Here we concentrate on the con- 
straints coming from the X-ray temperature and luminosity functions (XTF and XLF), 
and investigate several possible sources of discrepancies in the results. We conclude 
that the main source of error in both methods is the scaling relation involved, in 
particular the way its intrinsic scatter and systematic normalization are treated. 

For temperature derived constraints, we use a sample adapted from HIFLUGCS, 
and test for several sources of systematic error. We parameterize the mass-temperature 
relation with an overall factor T^,, which varies between about 1.5 and 1.9 in the 
literature, with simulations typically giving lower results than empirically derived 
estimates. After marginalising over this range of T*, we obtain a 68 per cent confidence 
range of erg = O-Tltool for a standard ACDM model. 

For luminosity derived constraints we use the XLF from the REFLEX survey 
and explore how sensitive the final results are on the details of the mass-luminosity, 
M-L, conversion. Assuming a uniform systematic uncertainty of ±20 per cent in 
the amplitude of the mass-luminosity relation by Reiprich & Bohringer, we derive 
as = 0.79t^;!5? ^r the same standard ACDM model. Although the XTF and XLF 
derived constraints agree very well with each other, we emphasize that such results 
can change by about 10-15 per cent, depending on how uncertainties in the L-T-M 
conversions are interpreted and included in the analysis. 

We point out that in order to achieve precision cosmology on erg using cluster 
abundance, it is first important to separate the uncertainty in the scaling relation 
into its intrinsic and overall normalization parts. Careful consideration of all sources 
of scatter is also important, as is the use of the most accurate formulae and full 
consideration of dependence on cosmology. A significant improvement will require the 
simultaneous determination of mass using a variety of distinct methods, such as X-ray 
observations, weak lensing, Sunyaev-Zel'dovich measurements and velocity dispersions 
of member galaxies, for a moderately large sample of clusters. 

Key Vifords: gravitation - galaxies: clusters: general - cosmology: theory - large-scale 
structure of Universe 



1 INTRODUCTION 

The normalization of density perturbations on large scales 
is a fundamental parameter describing our Universe. Its de- 
termination has been actively pursued for the last quarter 
century, using a wide range of methods. Now we are entering 
the era of precision cosmology, with several cosmological pa- 
rameters apparently determined to better than 10 per cent 
accuracy. While this level of precision has been reached for 
the amplitude of the largest scale perturbations, through the 



COBE anisotropy measurements, the direct normalization 
of matter fluctuations on galaxy cluster scales still yields a 
wide range of estimates. 

Here we focus on the determination of the cluster mass 
distribution via X-ray observations, and discuss how differ- 
ences in the details of the adopted strategies can lead to 
different, sometimes quite discrepant, determinations of erg. 
While considering X-ray observations, two approaches have 
traditionally been adopted. The first one is based on de- 
termining the cluster X-ray luminosity function (XLF) and 
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then converting it into the mass function through a suitable 
Mass-Luminosity, M-L, conversion (e.g. Henry et al. 1992; 
Sadat, Blanchard & Oukbir 1998; Reichart et al. 1999; Bor- 
gani ct al. 2001; Rcipricli & Bohringer 2002; Viana, Nichol 
& Liddle 2002; Schuecker et al. 2002). The second possibility 
is to use measurements of the X-ray temperature function 
(XTF) and then apply a Mass-Temperature, M-T, relation 
(e.g. Oukbir & Blanchard 1992; Eke et al. 1998; Viana & 
Liddle 1998; Pen 1998; Markevitch 1998; PierpaoU, Scott & 
White 2001; Seljak 2002). 

Measuring cluster luminosities requires ~ 10 times less 
X-ray photons than measuring temperature. This has al- 
lowed us to precisely determine the XLF for samples con- 
taining a few hundreds of local clusters. Current independent 
determinations of the XLF agree with each other at high pre- 
cision (see Rosati, Borgani & Norman 2002, for a review), 
thus making it a stable reference point to quantify the cluster 
population. However,the dependence of the bremsstrahlung 
emissivity on the square of the local gas density causes the 
X-ray luminosity to be highly sensitive to the core and to 
local details of the intra-cluster medium (ICM) structure. 
This makes it a noisy and difficult to calibrate mass estima- 
tor. Resorting to an observationally calibrated M-L relation 
is in principle possible through two alternative routes: ei- 
ther combining the L-T and the M-T relations, or directly 
measuring an M-L scaling. Both approaches require an in- 
dependent measure of M, which is a notoriously difficult 
problem. Traditionally, masses are estimated either by ap- 
plying hydrostatic equilibrium to a cluster with a measured 
temperature pattern and surface brightness profile (Reiprich 
& Bohringer 2002; Ettori, De Grandi & Molendi 2002), or 
using the mass estimated through cosmic shear (e.g. Viana 
& Liddle 2002, Allen et al. 2002). Both of these methods 
could suffer biases (Evrard, Metzler & Navarro 1996; Met- 
zler. White & Loken 2001) which need to be accounted for. 

X-ray temperatures are in principle easier to under- 
stand: since the ICM temperature is mainly determined by 
gravitational processes, it should be more directly related 
to the total gravitating mass of the cluster, i.e. the M-T 
relation should have less scatter and be more tractable the- 
oretically. Although this is true to a first approximation, 
the simplest expectation based on hydrostatic equilibrium 
and isothermal gas have been shown to provide a poor rep- 
resentation of the observed M-T relation (e.g. Finoguenov 
et al. 2001; Allen, Schmidt & Fabian 2001). In the absence 
of improved theoretical modeling one is forced to rely on 
an empirical calibration of the M-T relation, which raises 
many of the same issues as for the M-L relation discussed 
above. 

In general, all methods of mass estimation have both a 
significant scatter and potential systematic errors, leading 
to discrepant as results. 

For example, Pierpaoli, Scott & White (2001, here- 
after PSW) explored several improvements in determin- 
ing as from the distribution of local cluster temperatures, 
using a combination of X-ray data, cluster simulations 
and theoretical modeling. For = 0.3, they constrained 
the power spectrum normalization to lie around ag — 1, 
but stressed how existing uncertainties, particularly in the 
Mass-Temperature relation, dominate the error bar. 

Borgani et al. (2001; hereafter BOl, see also Rosati 
et al. 2002) analyzed the luminosity distribution for the 



ROSAT Deep Cluster Survey out to « ~ 1. They paid par- 
ticular attention to the impact of uncertainties and scatter 
in the M-L conversion on the final results. They preferred a 
lower normalization, as — 0.7 for Qm ~ 0.3, with a roughly 
15 per cent uncertainty mainly driven by the systematics 
affecting the M-L relation. 

In this paper wo use the local XLF from the REFLEX 
survey (Bohinger et al. 2002) and the distribution of tem- 
peratures from a new compilation of nearby clusters 

based on HIFLUGCS (Reiprich & Bohringer 2002). Be- 
sides providing further constraints on as, we will discuss how 
these two different approaches to the cluster mass function 
are afi'ected by systematics, and under which conditions they 
provide consistent constraints on cosmological parameters. 
Our basic conclusion will be that the two approaches in fact 
give consistent results, although at present it is not realistic 
to expect a determination of as to better than 10 per cent. 

More specifically, after giving the theoretical frame 
(Section 2), we analyze the temperature-derived constraint 
in Section 3. We assess several possible different sources of 
error, including: the M-T normalization and its scaling; the 
effect of the L-T normalization error on the computation 
of the effective volume; different temperature determina- 
tions; and different cuts of the data. For luminosity-derived 
constraints, in Section 4 we compare the results obtained 
from the local REFLEX sample with the RDCS, therefore 
testing possible evolutionary effects. We also compare re- 
sults obtained by applying the direct L-M relation with 
those obtained via the the L-T-M approach. In the end, we 
also compare the luminosity- and temperature-derived con- 
straints. We thoroughly discuss how error treatments can 
influence the final result throughout, particularly in Sec- 
tion 5, where we also compare our results to other recent 
determinations of as- Finally, we present our conclusions in 
Section 6. 



2 THEORY 

The mass variance is constrained from the cluster abundance 
through the mass function - the (comoving) number density 
of objects of a specified meiss. While there has been signifi- 
cant progress in the theory of the mass function of late, there 
still remain serious systematic uncertainties in the theoreti- 
cal predictions, and it is therefore important to specify which 
expression for the the mass function is being adopted. Jenk- 
ins et al. (2001) used a variety of iV-body simulations for 
different cosmological models to derive analytic expressions 
for the mass function of dark matter halos. In particular, 
the expression in their equation (B3) , 

rfM = °-3°^M^^fhrAr°^p(-i'^^" )' 

has been shown to reproduce the distribution of halos, if the 
mass is interpreted as that within the radius, risom, inte- 
rior to which the mean overdensity is p = 180pm, for their 
'tCDM' simulation. Here pm is the 'background' density or 
the mean cosmic mass density (i.e. f2m times the critical 
density). They also show that the mass function is approxi- 
mately 'universal' if the mass is taken to be the sum of the 
particles found in their simulation with a particular group 
finder. The combination of these findings strongly suggests 
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that equation (1) describes the mass function for a wide 
range of theories if the mass is measured interior to risom 
(see also Hu & Kravtsov 2002; White 2002). 

Jenkins et al. (2001) have confirmed and extended ear- 
her work which shows that the standard Press-Schechter 
expression for the mass function (Press & Schechter 1974) 
significantly overestimates and underestimates the number 
density of halos in the low- and high-mass end of their distri- 
bution, respectively (see also Evrard et al. 2002). They also 
found that the mass-function by Sheth & Tormen (1999) 
somewhat overestimates the number density of very massive 
halos, although this difference is probably negligible for our 
current purposes. However, in general, different expression 
for the mass functions induce differences of 4-8 per cent in 
the resulting value of ag , which are comparable to or larger 
than the statistical uncertainties. 

In addition to this, Jenkins et al. (2001) and White 
(2002) have shown that the mass function is not precisely 
a universal function, i.e. it is not simply a function of the 
linear density field smoothed on an appropriate scale. While 
different cosmologies predict mass functions which are very 
similar, in scaled units, there is a 'scatter' of several tens of 
per cent in number at fixed mass. The non-universality of the 
mass function introduces an additional, and non-negligible, 
systematic uncertainty in the conversion between a mear 
sured object's abundance and the cosmological parameters. 
As we push to more precise estimates of increasing num- 
bers of parameters these issues will need to be confronted 
head on. However, for the moment we shall ignore such re- 
finements as our focus will be on currently more important 
systematic errors. 

The mass entering in equation (1) should be interpreted 
as the mass contained within risom, i.e. inside a radius en- 
compassing a mean overdensity p = 180pm. However, scal- 
ing relations connecting mass to X-ray observable quantities 
may provide the mass at different values of p/p (see below). 
In this case we follow White (2001) and rescale the masses 
assuming an NFW profile for the DM halo (Navarro, Frenk 
& White 1996) with a concentration c = 5 appropriate for a 
rich cluster. A similar procedure by Hu & Kravtsov (2002) 
gives essentially identical results. 

The (juantity a{M,z) in equation (1) describes the per- 
turbation of the mass-scale M, which is the mass contained 
within a top-hat sphere of comoving radius R, such that 
M = {4-!r/3)pR^. In terms of the power spectrum, it is 

a\R, = j^^ ^'ik, W'ikR) , (2) 

where = P{k, z) /{2-k^), P(k) = |4P cc k"T^{k) is the 
matter power spectrum (n = 1 will be assumed in the follow- 
ing for the primordial spectral index), W{kR) is the window 
function corresponding to the smoothing of the density field 
(the Fourier transform of a top-hat in this case) and T{k) 
is the transfer function. For the latter we use the fitting 
expression provided by Eisenstein & Hu (1999). 

We notice here that some commonly used approxima- 
tions for as or for the transfer function may introduce sys- 
tematic sources of error. For example, the Viana & Liddle 
fitting formula for as as a function of cosmology is impre- 
cise at the level of 2-3 per cent for fim values as low as 
0.2. Another common choice is the use of the Bardeen et 
al. (1986) fitting formula with no baxyons in conjunction 



with an expression for the 'shape parameter' F which models 
the baryon dependence. This can introduce differences in the 
transfer function of 5-20 per cent around k ~ 0.1 /iMpc"^. 
We will stick with the Jenkins et al. (2001) mass function 
and the transfer function given by Eisenstein & Hu (1999). 

In carrying out these comparisons, we restrict ourselves 
to flat cosmological models with an initial power spectrum 
index of n = 1, and pure cosmological constant making up 
the dark energy. We will also fix the baryon fraction to be 
Ob/i^ = 0.02, and either allow fim to be free, or fix its value 
at 0.3. Scaling relations for other cosmologies can be found 
in PSW. 

Once one has adopted a cosmology and a mass func- 
tion, then the power spectrum can be normalized using some 
property of the clusters which can be used to estimate mass. 
For studies which focus on the X-ray properties of clusters, 
then all methods effectively use either the X-ray temperature 
function (XTF) or the X-ray luminosity function (XLF). 



3 CONSTRAINTS FROM THE XTF 

In this section we present our tjg determination derived from 
the temperature measurements. We first discuss the sam- 
ple used (Section 3.1), then the possible choices for tem- 
perature modeling (Section 3.2), the M-T scaling relation 
(Section 3.3), our methodology (Section 3.4), and finally the 
results (Section 3.5). 

3.1 The local sample 

Since almost all recently used cluster samples derive from 
essentially the same ROSAT databases, it is intriguing that 
different selections and treatment of the data apparently 
give such different results. It is therefore clearly important 
to define one's sample as carefully as possible. 

The error bars on as are not dominated by Poisson fluc- 
tuations in the small number of sources, but by systematic 
sources of scatter. So it is important to have a well under- 
stood cluster sample, rather than the largest available one. 
We speciflcally focus on a local sample, by which we mean 
clusters with red shift z < 0.1. This allows us to avoid evolu- 
tion issues. We also exclude the nearest clusters (at z < 0.03) 
to avoid possible biases introduced by survey incompleteness 
and local super-clustering. 

In PSW we deflned a cluster catalogue adapted from 
the approach of Markevitch (1998), which was based on 
-ROSAT-selected clusters with ASCA-dcrived temperatures. 
Our temperatures there were taken from White (2000), who 
fltted ASCA data with a model which included the effects of 
cooling flows. We supplemented these cooling-flow corrected 
temperatures with a few temperature estimates drawn from 
the literature. Since that work the available X-ray catalogues 
have improved significantly. 

The HIFLUGCS catalogue (Reiprich & Bohringer 2002) 
has now been published, and is probably the most complete 
publicly available X-ray catalogue based on the ROSAT 
All Sky Survey (RASS, Triimper et al. 1990). Reiprich & 
Bohringer (2002, hereafter RB02) discuss the selection of 
the HIFLUGCS sample in detail. Here we give the most rel- 
evant features only. The main selection criterion is a fiux 
threshold at 1.7 x lO'^Wm'^ in the ROSAT 0.1-2.4 keV 
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band. The final flux limit is set at /um = 2.0 x 10"" Wm"^, 
above which it appears the sample is very close to complete 
over the 8.14 sr it covers. This sample contains 63 clusters, 
with an additional 43 clusters in an extended list. Rather 
than use the ASCA derived, cooling flow corrected temper- 
atures from White (2000), we have used the temperatures 
quoted by Ikebe et al. (2002). They present temperatures for 
the HIFLUGCS clusters determined from a two-temperature 
flt. 

In constructing our sample we exclude both the near- 
est (a < 0.03) and most distant (z>0.1) clusters, to avoid 
possible biases. Within these limits wc wish to construct an 
effectively temperature limited sample, so wo need to adopt 
a relation between X-ray temperature and luminosity. The 
luminosity-temperature relation for HIFLUGCS is (Ikebe et 
al. 2002, equation (4)) 



L = (1.38 X 10^5) T^-s W, 



(3) 



with T measured in keV and L in the 0.1-2.4 keV energy 
band. 

Wc include a cluster if it would have passed the flux cut 
had its luminosity been given by L{T) and if it is hotter than 
3 kcV. The errors on the measured temperature are included 
by weighting each cluster based on the probability it is in the 
sample assuming a Gaussian distribution of temperatures 
described by the observed central temperature and error. 

While we use the mean L-T relation to compute the 
weight of each cluster, when we later compute the effective 
distance to which a cluster of temperature T could have been 
seen, we include the effect of scatter in L—T (which slightly 
reduces the distance). 

In the following we will use the primary HIFLUGCS 
sample (63 clusters) as the reference one, but we will also 
compare the results derived using the PSW catalog (72 clus- 
ters) and an extended version of the HIFLUGCS sample 
which contains 96 clusters from a combination of the ex- 
tended HIFLUGCS and PSW. Note that after applying the 
selection criteria presented above, the actual number of clus- 
ters used to determine as is significantly lower (see Sec- 
tion 3.5). 



3.2 Modeling cluster temperatures 

Authors differ in what they moan by cluster 'temperature' 
and there are a number of subjective choices which can be 
possible sources of discrepancy, including: different models 
for fitting T; varying fields of view and spatial resolutions; 
spectral band differences; and methods for dealing with sub- 
structure, modeling mergers, etc. The basic point is that the 
quoted errors do not include all of the systematic effects and 
care must be taken to avoid comparing apples with oranges. 

A particular difficulty is deciding how to deal with ra- 
dial structure, particularly in the cluster center. Clearly 
some procedure needs to be adopted, since a single isother- 
mal model does not adequately fit both the central and outer 
regions of many clusters. So long as the same procedure is 
adopted for the method used to fix the M-T normalization 
£is is used for the data, it probably makes little difference pre- 
cisely what is done. However this is rarely the case in prac- 
tice. In PSW, for example, we used temperatures from White 
(2000) which came from a specific cooling flow model which 
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Figure 1. Ikebe et al. (2002) isothermal temperature fits from 
ASCA data compared with de Grandi & Molendi (2001) temper- 
atures derived from Beppo-SAX data for those clusters in com- 
mon. This is indicative of the current size of systematic errors in 
temperature estimates. 



had been fitted to ASCA data, while Markevitch (1998) sim- 
ply excised the centres of clusters. The justification for this 
is that most simulations possess neither the resolution nor 
the physics to adequately model cooling fiows, and hence 
they should not be allowed to affect the data either. 

Recently Ikebe et al. (2002) performed a different fit- 
ting for the temperatures of the clusters. They use a 2- 
temperature picture, in which the outer regions are mod- 
eled by a single temperature isothermal plasma, while in 
the central regions there is a second cooler component form- 
ing a multi-phase plasma with the first. For some clusters, 
with no statistical measurement of a cooler component, they 
flxed its temperature at Thot/2. They argue that overall this 
method provides a reasonable description for both cooling- 
flow and non-cooling-flow clusters. Wc will therefore use the 
temperature of the hot isothermal gas as the mass proxy for 
each cluster in our sample. It is worth pointing out that the 
temperatures derived using this method axe typically lower 
than those obtained by some previous authors (see fig. 1 
in Ikebe et al. 2002); in particular, the discrepancy is big- 
ger for hotter clusters. For clusters hotter than 6keV, the 
White (2000) temperatures arc approximately 25 per cent 
hotter than those of Ikebe et al. (2002). It now seems likely 
that the White (2000) temperatures were overestimated in 
several cases. 

However, to further illustrate this issue we plot in Fig. 1 
a comparison of Ikebe et al. (2002) temperatures with those 
which are in common with a study of clusters using Beppo- 
SAX data by de Grandi & Molendi (2002). Firstly, it is clear 
that these entirely independent determinations of tempera- 
ture are in rather good agreement. However, in detail there 
are differences for individual clusters, and a general trend 
for the Ikebe et al. (2002) temperatures to be cooler for the 
highest T objects. This comparison should serve as a useful 
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guide to the size of possible systematic effects in the tem- 
perature determinations. We will later consider the effect on 
our as estimate of boosting the temperatures of the highest 
temperature clusters. 

In the following we will adopt temperatures from Ikebe 
et al. (2002) whenever available, that is for the main HI- 
FLUGCS sample, and the HIFLUGCS extended sample, 
while we adopt temperatures from PSW for the PSW sample 
and for the clusters in that supplementary sample. 

3.3 The Mass— Temperature relation 

The M-T relation is generally assumed to be: 

where T is in koV*, Ac is the mean ovcrdcnsity inside 
the virial radius in units of the critical density and = 
nm(l + 2)^ + nA-|-nk(l-|-2)^. Note that Ac depends on cos- 
mology, with fitting formulae given in PSW, and is a redshift 
dependent variable, which should be evaluated using the ap- 
propriate Q,k{z) and f2m(2). 

While equation (4) comes from rather simplistic argu- 
ments (essentially dimensional analysis and an assumption 
that clusters are self-similar) both observations and simula- 
tions suggest it is a good approximation for systems hotter 
than about 3 keV. For example Finoguenov et al. (2001) find 

Msoo = (4.22l°J^) X lO^^ft^o'M© x fcT^^'^io-l", (5) 

for kT > 3keV clusters, where Msoo is the mass within 
an overdensity of 500 with respect to critical. Other re- 
cent determinations are also consistent with a slope of 1.5 
(e.g. Allen et a. 2001; Xu ct al. 2001) provided that only the 
higher temperature clusters are considered. 

Unfortunately the level of agreement on the overall nor- 
malization, r, , is not so good (see e.g. Table 1 of Muanwong 
et al. 2002 or fig. 2 of Huterer & White 2002). In fact the 
uncertainty in T, has for some time been the dominant un- 
certainty in determining as- In PSW T« was assumed to be 
1.3, in BOl it was assumed to be 1.2. These values are rep- 
resentative of the older simulation results. The observations 
have almost always suggested a lower normalization for the 
M—T relation than implied by the theoretical model (i.e. a 
higher value of T*), which therefore leads to a lower as value. 
A rough scaling argument suggests (Huterer & White 2002, 
see also Fig. 2) 

ar^8 oc (T.)-°-« , (6) 

with a slightly weaker dependence if only the most massive 
clusters are used in the analysis (Evrard et al. 2002). In- 
cluding additional physics in the simulations and trying to 

* In PSW we called this normalization j3, following some of the 
simulation papers. However, this caused some confusion because a 
different subset of authors (motivated by the isothermal P model) 
have defined a /3 which is inversely proportional. In addition (3 is 
often used to denote the outer slope of the emission profile of 
clusters. Defining T* in this way is much less ambiguous. 



estimate the mass in the same way as the observers helps to 
reconcile the discrepant normalizations (e.g. Muanwong et 
al. 2002) and suggests that T* may be closer to 1.6 than the 
values assumed earlier. Some direct observational determi- 
nations (which estimate the mass in different ways) suggest 
an even higher value (Finoguenov et al. 2001, Xu et al. 2001). 

We would like to stress that there arc two kinds of un- 
certainty in the M-T relation: an intrinsic scatter in the 
temperatures of clusters of a given mass which arises due to 
differences in their formation history; and an overall uncer- 
tainty in the normalization of the M—T relation. The first 
sort of uncertainty is entirely statistical in character. It needs 
to be considered, because it both moves the central value 
for as and contributes to the error bar. Taking this scatter 
into account is straightforward if its distribution is known. 
Most studies have agreed that the scatter is around 10 15 
per cent in T, and is essentially a reflection of the different 
merging histories. The second sort of uncertainty is entirely 
systematic, however, and has to be treated in a different 
way, for example by marginalising the likelihood of ug over 
some prior distribution of values of T*. This distinction is 
not typically made in determinations of as from X-ray data. 

It is very easy for errors in the determination of T, to 
creep in. Recall that we are ultimately interested in deter- 
mining the mass within a sphere of radius rigom centered 
on the cluster in order to compare to theoretical mass func- 
tions. For a typical rich cluster risom ~ 2 — 3/i~^Mpc, well 
beyond the region where observations probe directly, with 
the possible exception of lensing studies. 

In order to determine this mass we may estimate Ma 

at one of various different density contrasts A using a va- 
riety of methods each known to have biases. Assuming an 
NFW profile with a constant concentration (c = 5 is as- 
sumed for cluster scale halos) and a flat cosmology with 
flm = 0.3, then some commonly used masses are related 
by Mvir = Ma, = 1.22 M200, and M500 = 0.72Af2oo, while 
Ac — 100. Hence wo can derive Mvir = 1.69 M500, so the 
M-T relation of equation (5) implies a normalization of 
T* ~ 1.9. Significant corrections arc being applied in these 
conversions, involving large extrapolations from given pro- 
files, indicating that there is plenty of room for error. 

We caution that there is also a significant scatter in 
cluster concentrations, and that the NFW profile is only 

'universal' in that it fits statistically to isolated clusters, 
rather than holding in detail for each cluster. We believe 
that these issues do not have a large effect on determining 
the cluster normalization (White 2002), however this could 
be studied in more detail. Careful direct comparison between 
simulations and observations is probably ultimately the best 
approach when the next generation of simulations becomes 
available. 

In summary, while both observations and simulations 
agree quite well on the scaling law of equation (4) , there is no 
obvious agreement on the value of the normalization T* . The 
observations may favour a higher T«, perhaps 1.8-1.9, while 
the simulations at present appear to give values around 1.5- 
1.6. At present it is by no means clear which should be the 
preferred value, or indeed if they are genuinely discrepant. 
In what follows, we will investigate the implications on the 
determination of as of adopting different M-T assumptions 
and error treatments. 
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3.4 Methodology 

The general method is similar to the one described in PSW, 
with some technical improvements. We developed two in- 
dependent codes for the likelihood evaluation, with slightly 
different approaches for the treatment of the errors, and wc 
cross-checked the results. This allows us to to identify areas 
of potential disagreement in the erg estimate. 

Perhaps the largest modification is that our former 
Monte Carlo method for the error treatment has been com- 
pletely replaced by 'weight functions' to increase speed. For 
each set of cosmological parameters, including the normal- 
ization and slope of the mass temperature relation, we con- 
struct a grid of n{T)dT, taking into account the asymmetric 
nature of the quoted observational errors on the tempera- 
tures. At this stage a volume cut based on the L T relation 
is applied. A conversion from temperature to mass, and be- 
tween mass definitions, is used to compare the observed n(T) 
with theory using the Poisson averaged likelihood 

{log £) = ^ di log Hi- tJbi, (7) 

where di is the fraction of the cluster in bin i and jii is the 
mean value computed from the mass function and effective 
volume. 

In principle the uncertainty in the mass temperature 
relation could be treated in two different ways: either with 
a convolution of the mass function or by marginalising the 
value of the likelihood. The former approach is appropri- 
ate for treating intrinsic scatter, while the latter is more 
appropriate for an uncertain parameter such as the overall 
normalization or slope. In PSW wc showed that the un- 
certainty in relating the mass to the observable dominates 
the final error bar, so it is important to treat this issue as 
carefully as possible. Both approaches have been adopted in 
the literature, accordingly, we investigated the influence of 
the two approaches on the likelihood independently. Wc also 
checked that the intrinsic scatter could be treated equiva- 
lently through either convolution or by an extended Monte 
Carlo procedure. 

3.5 Results 

We present here our results, quoted with respect to a ref- 
erence model which assumes the HIFLUGCS extended cat- 
alogue with temperature from Ikebe et al. (2002). Unless 
otherwise stated, we will focus on T* = 1.75 and fim/i = 0.2. 
Our fiducial choice for scatter in the M~T relation is 10 per 
cent in T, with a Gaussian distribution. 

We show in Fig. 2 the variation of ag with T, if wc 
hold fim = 0.3. We find that the uncertainty due to the 
choice of the fiducial T» value dominates with respect to 
the assumed intrinsic scatter. In order for this situation to 
be improved in future studies, the discrepancy between the 
observational and simulated T* values will need to come into 
better agreement. 

In Fig. 3 we show the effect of treating the uncertainty 
in T, in two different ways: via convolution of the likelihood 
function or of the mass function. It is apparent that after 
marginalisation of the likelihood is performed, the overall 
effect is to broaden the error bars with respect to the case 
with no errors. Convolution of the mass function, instead. 



mainly shifts as toward lower values. The two approaches 
are therefore not equivalent, and should be used appropri- 
ately (i.e. convolution or Monte Carlo for the intrinsic scat- 
ter and marginalisation for the systematic normalization un- 
certainty) . 

These results have been obtained with Gaussian scatter 
in T. Some authors (e.g. BOl) have considered the scatter 
instead to be Gaussian in log(T), as observed scaling rela- 
tions sometimes quote. While there is no particular reason 
for the latter choice, apart from the convenience of fitting 
a straight line to the data, we verified that the two differ- 
ent choices give similar error bars, at least at the la level, 
and for the 10 per cent scatter assumed here. However, we 
caution that in practice the distribution may not be Gaus- 
sian, and the tails may have a disproportionate effect. We 
also note in Fig. 3 that for relatively small intrinsic error (10 
per cent in T) the effect of the scatter is sub-dominant with 
respect to the other sources of error. However, we see that 
a 16 per cent scatter is sufficient to shift the as value by a 
significant amount. 

We show in Fig. 4 the dependence on the assumed mass 
function. We note here that the results derived with the 
Jenkins at al. (2001) mass function are more similar to the 
ones obtained using Press-Schechter (1974) than to those 
assuming the Shetli & Tormen (1999) form. This is due to 
the fact that the Jenkins et al. (2001) function is actually 
closer to Press-Schechter than to Sheth-Tormen in the range 
0.3 < lna~^ < 0.7 which is mostly relevant here. 

The dependence on the specific sample used is shown in 
Fig. 5. We note that although the cluster sample contains 
97 clusters in its extended version, due to our selection cri- 
teria the actual number of clusters used is much lower: the 
code selects on average 28.5 clusters out of the HIFLUGCS 
sample, 38 out of the extended one and 41.5 out of the old 
PSW one. About 10 per cent discrepancy between the HI- 
FLUGCS sample and the PSW one is due to the inclusion 
of the clusters that arc nominally below the flux limit (and 
scatter up from the extended sample). These clusters there- 
fore give a non-negligible contribution to the final as value 
through our likelihood technique. The rest of the discrep- 
ancy may be explained with the different modeling of the 
temperatures between Ikebe et al. (2002) and White (2000). 

For the dependence of as on cluster temperature, we 
verified that cutting out clusters below 4keV does not re- 
sult in a significant difference. This means that our 3keV 
low temperature cut is high enough to include clusters that 
actually follow the same M-T relation. We also boosted the 
temperatures of these clusters above 6.5 keV by 10 per cent, 
finding no noticeable change in the 90 per cent contour lev- 
els. This implies that the level of discrepancy in temperature 
estimation between, for example, Ikebe et al. (2002) and de 
Grandi et al. (2002) is probably not a major source of un- 
certainty. 

We also tested the error in the slope of the M-T re- 
lation, finding that it is not very important. Increasing the 
exponent from 1.5 to 1.6 lowers as by 0.05, while decreasing 
it to 1.4 increases as by a similar amount. Since we define 
our pivot point at 10^^ h~^yiQ, above where most of the 
sample lies, the sense of the shift is as expected. 

Finally, introducing a 35 per cent scatter in the L-T re- 
lation when computing volume implies a reduction of about 
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Figure 2. The dependence of trg on the value of the M-T nor- 
malization T» assumed. Here we fix rJni = 0.3 and h = 0.7. Error 
bars are la assuming an intrinsic scatter of 10 per cent in T» . The 
choice of T* is probably the most important difference between 
some discrepant values of as in the literature. 



Figure 4. The 90 per cent confidence contours obtained with 
different mass functions for a fixed T, = 1.75. The solid line 
corresponds to the Jenkins et al. (2001) mass function in equa- 
tion (1), the dashed line to that of Sheth & Tormen (1999), and 
the dot-dashed to the Press & Schechter (1974) formula. 




O. 60 



Figure 3. 90 per cent confidence regions in the (trs.rJi^i) plane 
from the analysis of cluster temperatures. The solid contour cor- 
respond to a fixed T, = 1.75 (no T scatter), the dot-dashed is 
obtained by convolving the likelihood function with a 10 per cent 
error in T, while the three-dots-dashed and dashed are obtained 
by convolving the mass function with a Gaussian window with a 
10 per cent and 16 per cent width in T. 



Figure 5. The 90 per cent confidence contours obtained using 
different input catalogues for a fixed T« = 1.75. The solid line 
is obtained using the HIFLUGCS sample, the dot-dashed with 
the extended version (the most appropriate sample to use in our 
particular methodology), and the dashed with the old PSW cat- 
alogue. In the latter case, the L-T relation used is the one from 
PSW. 
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8 per cent in the range 4 < T < 8 keV. However, this has a 
neghgible effect on the final as result. 

After investigating all of these effects, it is clear that 
one has to be very careful in order to obtain a robust es- 
timate for (78 • But even when that has been achieved, the 
overall error is dominated by the systematic uncertainty in 
the normalization of the M-T relation. 

To obtain an overall 'best efforts' estimate for we 
can marginalise over a uniform distribution for T, , stretch- 
ing between typical simulation values (1.5) to typical obser- 
vational values (1.9). For the standard ACDM model (with 

= 0.3, h = 0.7, Q.bh^ = 0.02 and n = 1) we obtain 

= 0.771°;°! (8) 

(c.f. Fig. 2), where the error bars describe the 68 per cent 
confidence region (corresponding approximately to ila). 
Here wo use our best attempt to construct an effectively 
temperature-selected sample, which consists of the HI- 
FLUGCS extended sample (with suitable cuts), plus some 
additional clusters from PSW. Scalings similar to those given 
in PSW could be used to extend this to different cosmolo- 
gies. 



suits have been also found by Allen et al. (2002), who used 
the XLF from the REFLEX survey and the BCS (Ebel- 
ing et al. 2000), combined with an L-M relation calibrated 
from Chandra/ ROS AT data and weak lensing observations. 
Rather than focusing on a further determination of cosmo- 
logical constraints from the REFLEX XLF, our analysis is 
more aimed at understanding how such constraints are af- 
fected by the details of the analysis method, and how they 
compare with those derived from the XTF analysis, pre- 
sented in the previous section. 

Since we use the binned representation of the XLF given 
by B02, wo determine best-fitting values for cosmological pa- 
rameters and confidence levels by computing the between 
the measured and predicted REFLEX XLF. The model XLF 
is computed from the mass function of equation (1) by con- 
verting masses into observed X-ray luminosities (see be- 
low). After finding the values of the {ag,QYn) parameters 
which minimize the value, confidence regions arc esti- 
mated from standard 2-dimensional Ax^ variations. We as- 
sume the REFLEX XLF to be determined at the effective 
redshift Zgs = 0.1, which is close to the median redshift of 
the survey. 



4 CONSTRAINTS FROM THE LOCAL XLF 

An alternative approach is to dispense with the X-ray tem- 
peratures and to use the X-ray luminosity function (XLF). 
This method suffers from the potential limitation that the 
relation between X-ray luminosity and mass is less well un- 
derstood and is more affected by uncertainties in the ther- 
modynamic status of the ICM than the M-T relation. How- 
ever, the XLF has the significant advantage that it is now 
precisely calibrated for extended sets of both nearby and 
distant {z ~ 1) cluster samples. Quite remarkably, all the 
determinations of the local XLF based on the ROSAT All 
Sky Survey (RASS) and on deep pointings agree pretty well 
with each other(see Rosati, Borgani & Norman 2002, for a 
review), thus providing a stable reference point to establish 
the evolution of the X-ray cluster population. 

Therefore, it is important to understand the extent to 
which XTF and XLF approaches give consistent results on 
cosmological constraints. This is only possible through a 
careful consideration of the different steps involved. The 
most recent determination of the local XLF is based on the 
REFLEX survey (Bohringer et al. 2002, B02 hereafter). This 
survey includes about 450 clusters out to z <^ 0.3, and is 
complete down to the fiux limit of 3 x 10~^^Wm~^ in the 
0.1-2.4 keV energy band. B02 provided a binned representa- 
tion of the cluster XLF in the luminosity range L ~ 4 x 10^^ 
to 3 X lO^'^W (with h = 0.5) for both an Einstein-de Sitter 
cosmology and for a flat low-density model with Qm = 0.3. 
The effect of changing cosmology is almost negligible at the 
low redshifts probed by the REFLEX survey. For definite- 
ness, we will use in the following the XLF determination for 
the low-density model. 

Schuecker et al. (2002) combined clustering properties 
and redshift distribution of REFLEX clusters, through a 
mciximum-likelihood approach, to obtain constraints on as 
and fim. They applied the L-M relation calibrated by RB02 
and obtained its — 0.75 for fim = 0.3 with formal sta- 
tistical uncertainties of about 5 per cent. Consistent re- 



4.1 Converting masses into X-ray luminosities 

In order to convert the observed L into mass to use in the 
mass function of equation (1), we follow two different pro- 
cedures. 



4.1.1 Method (a): using L-T and M-T 

The first method combines the observed L-T relation and 
the M-T relation calibrated from hydrodynamical cluster 
simulations. This procedure has been applied by several au- 
thors in the analysis of flux-limited cluster surveys (e.g. Ki- 
tayama & Suto 1998, Sadat, Blanchard & Oukbir 1998, Re- 
ichart et al. 1999, BOl). As for the M-T relation, we use the 
expression of equation (4). 

As for the L-T relation, different analyses consistently 
show that is well represented at cluster scales (T ^ 2 keV) by 
a power law, L ocT", with a ~ 3, and with a scatter which is 
significantly reduced once the effect of cooling flows has been 
corrected for (e.g. Markevitch 1998, Allen & Fabian 1998) or 
by considering only systems without significant cooling-flow 
signatures (e.g. Arnaud & Evrard 1999). Following the same 
notation for the relation between bolometric luminosity and 
emission- weighted temperature as in BOl, the expression we 
use is 

Lboi = L6(r/6keV)" W^\~'^ W , (9) 

with L(i = 3 and a = 3. Bolometric luminosity is then con- 
verted into the 0.1-2.4 keV band by using a mekal spectral 
synthesis code, assuming one-third solar metallicity for the 
intra^cluster medium. Since this conversion has a larger ef- 
fect for hotter systems, the L-T relation tends to flatten 
when computed in the 0.1-2.4 keV band. This is the rea- 
son for the different slopes in equations (3) and (9). We 
assume that the combination of M-T and L-T variation re- 
sults in an overall 45 per cent scatter, which is assumed to 
be Gaussian-distributed in the log. 
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4.1.2 Method (h): using M-L directly 

The second method uses an observationally determined 
M-L relation (RB02, Ettori et aL 2002). RB02 estimated 
the M-L relation for a sample of 103 clusters with L ^ 
10^^/iJg^W, temperature measured from ASCA and gas den- 
sity profiles from ROSAT PSPC pointed observations. Et- 
tori et al. (2002) used Beppo-SAX data for a set of 22 clus- 
ters with T > 3 keV, with measured profiles of gas tem- 
perature and density. Based on such data, they found that, 
unlike for the L-T relation, no segregation between cooling- 
flow and non cooling-flow clusters exists in the M-L plane. 
As a consequence, these authors claim that directly mea- 
suring the M-L relation is more stable than combining the 
L-T and M-T relations. 

In the following, we will use the M-L relation from 
RB02, which includes lower-luminosity systems, thus bet- 
ter covering the L range probed by the REFLEX XLF. 
We consider the relation based on mass estimated at over- 
density A — 500 and luminosity in the 0.1-2.4 band. As 
in Section 3.3, M500 will be related to Ma^ assuming an 
NFW profile with c = 5 for the concentration parame- 
ter. RB02 claims that no extrapolation is required for 86 
per cent of the clusters in their sample to obtain mass 
at such an overdensity (although obviously the assumption 
of a c = 5 NFW profile is still required to connect M500 
to Misoni). RB02 adopted a fitting expression of the form 
L = 10-^ (M5oo//i^o^MQ)"/ig[flO^^W and provided fitting 
parameters for both the direct {A = —19.708, a — 1.652) 
and inverse relations {A = —17.545, a = 1.504; see Table 
10 in RB02). In the following we will show results based 
on these different choices for A and a, so as to judge the 
stability of constraints on cosmological parameters against 
uncertainties in the fitting procedure. 

Finally, RB02 also estimated the overall scatter in mass 
from this relation to be ~ 48 per cent. After subtracting 
in quadrature the contribution from the mean mass mea- 
surement uncertainty, ~ 28 per cent, the resulting intrinsic 
scatter in the M-L relation turns out to be 39 per cent. 
This value overestimates the true scatter if errors in the 
meiss measurements have been underestimated, for instance 
due to the assumptions of isothermal gas and /3-model fit- 
ting for the gas density profile, made by RB02 in their anal- 
ysis. Also, if significant temperature gradients are present 
(e.g. Markevitch et al. 1998, Finoguenov et al. 2001, De 
Grandi & Molendi 2002, Pratt & Arnaud 2002), then the 
isothermal assumption would lead to an overestimate of the 
total gravitating mass and, therefore, to an overestimate of 
as at fixed f2m. Whether violations of these assumptions 
would lead to a bias or to an increase in the scatter or both 
is unclear. 

Following the same line of reasoning as in Section 3.4 for 
the M-T relation, any uncertainty in the M-L relation can 
be interpreted as due either to a genuine intrinsic scatter or 
to an overall uncertainty in the normalization of the relation. 
In the first case, one has to convolve the model mass function 
with a scatter, while in the second case the correct procedure 
is to marginalise over the amplitude of the M-L relation, 
after assuming a range of variation and a distribution for this 
parameter. In the following we will assume as our fiducial 
analysis that based on the RBQ2 M-L relation, with 20 per 
cent intrinsic scatter in mass. We will also verify how the 
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Figure 6. The relation between M500 and luminosity in the 0.1- 
2.4 keV energy band. Data points are from Reiprich &; Bohringer 
(2002, RB02). Solid and dotted lines represent the direct and 
inverse best- fitting relations from RB02. Long-dashed and short- 
dashed lines have been obtained from the combined L-T and 
T-M method to convert mass into L, after assuming T* = 1.75 
and T« = 1.2, respectively. We assume Hm = 0.3 and use equa- 
tion (4) for the M-T conversion. Bolometric luminosities have 
been converted into L in the 0.1— 2.4keV band by using a MEKAL 
model with Z = 0.3 Zq for the ICM metallicity. 

results change by marginalising over the amplitude of the 
AI-L relation and by changing the intrinsic scatter between 
the two extreme assumptions of negligible scatter and 39 
per cent scatter. We show in Fig. 6 the observational data 
points from RB02 along with the different M-L relations 
that we use in our analysis. Method (a) with T* — 1.2 (used 
by BOl) overestimates the mass at a fixed luminosity, while 
taking T, — 1.75 is rather close to the fitting relations used 
by RB02 for L ;> lO^^W. 

4.2 Results 

^From their analysis of the ROSAT Deep Cluster Survey 
(RDCS), BOl derived constraints on the (ilnncrs) plane by 
following the evolution of the cluster population out to z ^ 
0.8 (see also Rosati et al. 2002). This analysis was based 
on the combined L-T and T-M approach of method (a) 
with T, — 1.2, no significant evolution of the L-T relation 
(e.g. BOl, Holden et al. 2002, Novicki, Sornig & Henry 2002; 
cf. also Vikhlinin et al. 2002), an overall 45 per cent scatter 
in the M-L relation, and assuming the mass function by 
Sheth & Tormen (1999). 

For the sake of direct comparison, we repeat here the 
same analysis on the REFLEX XLF. We emphasize that a 
major difference between the analysis of RDCS and of RE- 
FLEX lies in the widely different volumes covered by the two 
surveys at different redshifts - unlike REFLEX, the RDCS 
has a modest volume coverage at low redshift, while being 
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Figure 7. Comparison between the constraints derived from the 
REFLEX and the RDCS XLF. The solid contour indicates the 90 
per cent confidence region from the REFLEX XLF, by using the 
method (a) for the mass-luminosity conversion (Section 4.1.1), 
with L-T relation given by equation (9), M—T relation by equa- 
tion (4) with T, = 1.2. The dotted contour corresponds to the 
analysis of the evolution of the RDCS XLF out to 2 ~ 0.8 (BOl, 
Rosati et al. 2002; see text), for the same choice of L—T and M- 
T relations and assuming a non-evolving L-T (for Qm = 1.0). 
In both analyses the mass function is convolved with 45 per cent 
scatter in the overall M-L conversion. For the sake of comparison 
with the results of BOl, we used here a BBKS power spectrum 
with r = 0.2 for the shape parameter and the mass function by 
Sheth & Tormen (1999). 



able to trace the cluster population out to « ~ 1. Therefore, 
the analysis of the RDCS provides a dynamical constraint 
on Qm, by probing the growth rate of density perturbations, 
while REFLEX provides a geometrical constraint on Qm 
through the shape of the power spectrum (assuming CDM 
and fixing h and fib from independent observations). 

We compare in Fig. 7 the constraints from REFLEX 
and RDCS. Quite remarkably such constraints are in good 
agreement. The somewhat looser constraints from REFLEX 
should be understood in terms of the largely different nature 
of the two samples, which cover widely different rcdshift 
ranges, and of the different analysis method - the RDCS 
analysis is based on a maximum likelihood approach, which 
allows one to extract the whole information provided by the 
distribution of clusters in the L-z plane. Besides the values 
of (Tg, the agreement between the 'geometrical' and 'dynam- 
ical' measurements of ilm is quite remarkable. We believe 
that this finding represents important support for the relia- 
bility of using galax:y clusters as tracers of the evolution of 
cosmic structures. 

We show in Fig. 8 how cosmological constraints alter as 
we change different aspects of the analysis procedure, such as 
the M-L conversion, as well as the amount of uncertainty in 
this conversion and how it is treated. Our reference analysis 
corresponds to choosing the best-fitting direct M-L relation 



from RB02, with 20 per cent intrinsic scatter, convolved with 
the mass function (results shown with solid contours in the 
panels). In this case, we obtain 



<78 
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(10) 



where error bars correspond to 68 per cent uncertainties 
on the two parameters. If we fix the density parameter at 

f2ni = 0.3, we find 



as = 0.74: 



(11) 



where the error corresponds to the 68 per cent uncertainty 
for one interesting parameter. The above quoted errors only 
reflect the shot noise, propagated from the Poissonian er- 
rors quoted by B02 (and no systematic uncertainty in the 
M-L relation, which we discuss below). As pointed out by 
Evrard et al. (2002; see also Hu & Kravtsov 2002 and White 
2002), sample variance due to finite survey size and cluster- 
ing should also be considered. However, given the large size 
of the REFLEX survey, we expect sample variance not to 
be important, compared with other uncertainties related to 
the M-L conversion. 

In the left panel of Fig. 8 we show the effect of chang- 
ing the amount of intrinsic scatter in the M-L relation. As 
expected, assuming a larger scatter decreases the normaliza- 
tion of the power spectrum. For instance, taking Qm ~ 0.3 
and assuming zero scatter gives as ~ 0.8, while assuming 
a 39 per cent scatter lowers as to about 0.65. In the cen- 
tral panel we show the effect of the statistical uncertainty 
in the M-L relation, by considering both the direct and the 
inverse fitting provided by RB02, as well as the effect of 
marginalising the likelihood over the M-L amplitude. The 
effect of using the inverse fitting (dotted contour) is that of 
marginally favouring lower fim values, while keeping the ag- 
Qm degeneracy direction unchanged. This indicates that the 
statistical stability of the M-L conversion has only a small 
effect. As for the marginalisation procedure (thin continuous 
contour), it has been realized by integrating the likelihood 
function over a range of M-L amplitudes, defined by vary- 
ing the best-fitting value by ±20 per cent, with a uniform 
distribution. The effect of marginalising is that of widening 
the contours, while shifting the most likely as to larger val- 
ues, since no convolution with the scatter is now performed. 
In fact, for fim = 0.3 we have 



as — U.<'9_o.o7 ■ 



(12) 



for the 68 per cent confidence region. 

The right panel shows the effect of using method (a) for 
the M-L conversion, and assuming T« = 1.75, as a compro- 
mise between observational and theoretical determinations. 
Since the M-L relations from these two methods are quite 
similar (see Fig. 6), the slightly higher normalization from 
method (b) is precisely the consequence of having assumed 
a 20 per cent intrinsic scatter, instead of 45 per cent as in 
method (a). 



5 DISCUSSION 

5.1 Comparing XTF and XLF results 

We can now compare the XLF and XTF approaches, mak- 
ing as far as is possible consistent assumptions throughout. 
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Figure 8. The 90 per cent confidence contours for the REFLEX XLF, for different clioices of tfic M-L conversion. In all panels, the 
solid contour corresponds to the reference analysis method, taking the best-fitting direct M500— i conversion by Reiprich & Bohringer 
(2002, RB02) and convolving the model mass function with 20 per cent intrinsic scatter. Left panel; the dotted contour is for the effect 
of assuming no intrinsic scatter in the M500-L relation by RB02, while the thin continuous contour is for the maximal assumption of 39 
per cent intrinsic scatter in the Msoo~L relation. Central panel: the thin continuous contour shows the effect of marginalising over the 
amplitude of the M-L relation, which is assumed to vary with a uniform distribution of ±20 per cent about the best-fitting value; the 
dotted contour is based on taking the inverse M-L best-fitting parameters, as provided by RB02, in the reference analysis. Right panel: 
the dotted contour is for the M-L conversion based on the same method used to obtain the solid contour in Fig.7, but with T* = 1.75 
for the M-T relation (long-dashed line in Fig. 6). 
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Figure 9. Comparison between the 90 per cent confidence con- 
tours derived from the REFLEX XLF (continuous line) and the 
XTF analysis (dotted line) . The XLF analysis is based on method 
(6) for the M-L conversion, by marginalising over the amplitude 
of the M-L relation by RB02, which is varied by ±20 per cent 
around the best fitting value with uniform distribution (see text). 
The XTF contour has been obtained by marginalising over T*, 
which is assumed to vary with uniform distribution from 1.5 
to 1.9 (see also eq.lO). Both analyses are based on assuming 
Ui, = 0.019 h^'^ and h = 0.7 in the power spectrum by Eisen- 
stein & Hu (1999). 



We show in Fig. 9 a comparison between the marginaUsed 
constraints derived from the analyses of the previous sec- 
tions. Both analyses have been realized by assuming fit = 
0.019 ft^^ and h = 0.7 in the power spectrum by Eisen- 
stein & Hu (1999). The XLF constraint is that derived by 
marginalising over the amplitude of the M-L relation by 
RB02, on which the resuh of eq.(15) is based. The XTF 
result is derived by marginalising over T, , which is uni- 
formly varied from 1.5 to 1.9. Even bearing in mind the 
systematics affecting both the XTF and the XLF analyses, 
the two methods provide remarkably consistent constraints. 
This represents a non-trivial result, owing to the quite differ- 
ent approaches by which such constraints have been derived: 
not only are the cluster samples and the analysis meth- 
ods different, but also the systematics associated with the 
M-T-L scalings enter in completely independent ways. Al- 
though different assumptions on how to interpret and treat 
the scatter in the M-T and M-L relations produce different 
XLF- and XTF-based constraints, nevertheless it is reassur- 
ing that such a good agreement is just obtained for rather 
natural choices of these parameters. The wider contour from 
the XTF analysis is due to the narrower dynamical range 
probed in cluster masses by our sample of cluster temper- 
atures, which does not allow to put strong constraints on 
the shape of the power spectrum. In fact, while our XTF 
samples about half a decade in temperature, the REFLEX 
XLF samples about 3.5 decades in cluster luminosity (0.1- 
2.4 keV band), which correspond to almost 1.5 decades in 
temperature from the L-T relation of eq.(3). 

The main disagreement between the XTF analysis by 
PSW and the XLF analysis by BOl was neither due to the 
reliability of one of the two or both samples, nor to the T, 
assumed (which in fact would have implied a discrepancy in 
the opposite direction), but rather to the procedure of con- 
volving versus marginalising the scatter. This illustrates that 
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although the appUcation of the cluster abundance normal- 
ization technique appears to be simple, it involves a chain 
of steps, and the source of discrepancies can be rather sub- 
tle. However it also serves to reinforce our central point - 
with sufficient care the statistical and methodological differ- 
ences can be reduced to the point that the major remaining 
uncertainty is in the normalization of the mass scale of clus- 
ters and in the determination of the nature of the scatter 
between the theory-predicted mass and the observed X-ray 
luminosity and temperature. 

5.2 Comparing with other estimates of as 

We list in Table 1 other recent determinations of as, which 
are based on different methods. Older estimates and associ- 
ated references can be found in PSW. Rather than analyse 
individual studies in detail, we will make some general re- 
marks which help explain the differences. 

The normalization of density perturbations (through 
the mass variance at 8 /i^^Mpc) can be evaluated by means 
of several different observables, either on much larger scales 
(e.g. CMB derived constraints), on scales ~ 10 h~^Mpc (con- 
taining the mass of a rich cluster of galaxies, which is not 
far into the non-linear regime), or on much smaller scales 
(e.g. through lensing or Ly a forest variance measurements). 
Normalization estimates from large or small scales can be in- 
terpreted as as determinations through extrapolation of the 
power spectrum. This is cosmology dependent, and hence 
it should be no surprise that some of the estimates in Ta- 
ble 1 appear discrepant. Even if we focus on methods which 
are directly measuring the mass variance on scales close to 
8/i~^Mpc, there still appears to be a bewilderingly wide 
range of published values. 

One point to realize is that different cosmological as- 
sumptions are made by different sets of authors. Since the 
density parameter is constrained to a first approximation 
through the combination fim/i (sometimes referred to as the 
'shape parameter'), amy uncertainty in the value of h trans- 
lates into a similar uncertainty when constraining Clm ■ Hence 
there are ambiguities here, which one has to be careful about 
when comparing results between authors. Different values of 
as will be obtained using identical data but either a fixed 
h, fixing Qmh, or marginalising over various priors. This ex- 
plains some of the discrepancies in the Table. 

There have been several recent results from weak lens- 
ing (Bacon et al. 2002, Brown et al. 2002, Hoekstra et 
al. 2002, Jarvis et al. 2002, Refregier et al. 2002, van Waer- 
beke et al. 2002). The quoted as values, reported in Table 
1 in chronological order, range between 0.7 and 0.9, and 
have been decreasing in the past year. It is early days for 
this method, since the formal significance of cosmic shear 
detection is still rather modest. More importantly for the 
purposes of comparison with cluster abundance values of 
178 , most lensing studies have only poor knowledge of their 
source redshift distributions and are measuring on a scale 
which is a factor of several smaller than the ~10/i~^Mpc 
probed by clusters, and may be affected by non-linear clus- 
tering. Despite the fact that latest works have done major 
improvements in assessing systematic error via the B-mode 
power spectrum, this still remains a very hard measurement 
and the nature of the systematic errors is not yet well under- 
stood. Estimates using large-scale structure measurements 



coupled with Cosmic Microwave anisotropics (e.g. Lahav et 
al. 2001) also give similar values, but are even more depen- 
dent on the adopted cosmology and which parameters are 
held fixed. 

There are many other minor differences between some 
of the studies involving the analysis of the cluster number 
density. Precisely how the samples are constructed is obvi- 
ously a major issue. For instance, when dealing with X-ray 
properties this is particularly an issue for deciding how best 
to deal with substructure and cooling fiows. Here different 
authors make different decisions, and it is to be hoped that 
such ambiguities will eventually disappear as the simulations 
improve and can be more directly compared with the obser- 
vations. Another issue is merging corrections, which some 
recent studies have used, motivated by the idea that Press- 
Schechter describes the number density of halos at the viri- 
alisation epoch. However (as we stressed in PSW) , it is clear 
that the effect of merger history is already implicit in the 
intrinsic scatter of the M T relation, and so should not be 
included twice. Other minor issues include whether redshift 
corrections are made, what profiles are used to convert be- 
tween different mass quantities, whether there is curvature 
in the adopted M-T or M-L relations, etc. 

To summarize this part: differences in as can be due to 
a number of factors. Firstly, the method used may require 
an extrapolation to 8/i~^Mpc. Secondly, there are varying 
assumptions made about the values of other cosmological 
parameters in the analysis. Thirdly, there are different forms 
used for the mass function or the transfer function. Fourthly, 
even when the basic method adopted is the same, there can 
still be detailed differences in the approaches of each set 
of authors (e.g. which sources of scatter are included, how 
the catalogue is constructed or whether a merger correc- 
tion is applied). And finally, for XTF- and XLF-based ap- 
proaches there is the choice for the mass-temperature and 
mass-luminosity conversion, respectively. 



6 CONCLUSIONS 

We investigated the determination of as as derived from the 
XTF and XLF approaches. We considered many possible 
choices in terms of data and methodology, and paid special 
attention to the assessment of possible sources of biases and 
systematic uncertainties in the adopted procedure. 

In both XTF- and XLF-derived constraints, the most 
important effect on as estimates is the error in the scaling 
relation connecting mass to temperature and to luminosity. 
Moreover, different ways of introducing such an error have 
a significant impact on the final as range. If it is introduced 
as a convolution with the mass function, its overall effect is 
to lower the as value, while if it is introduced by marginal- 
ising the final likelihood, then the consequence is mainly to 
broaden the error contours. We argue that the former (lat- 
ter) procedure is more adequate when dealing with statisti- 
cal (systematic) errors. For instance, in the marginalisation 
procedure one has to make a decision about the functional 
form of the distribution with which the marginalised pa- 
rameters are varied. In our comparison between XTF- and 
XLF-based constraints (see Fig. 9) we decided to assume 
a uniform distribution for T* within a finite range. While 
there is no rigorous justification for this choice, this pro- 
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Table 1. Some recent estimates of erg. The errors are typically statistical. We have tried to give values consistently for — 0.3, while 
r ~ 0.2 is a typical value of the shape parameter. Methods used include: weak lensing (WL); Optical clusters (OC); clusters normalized 
using weak lensing (WLC); X-ray temperature function (XTF); power spectrum (PS); Sunyaev-Zel'dovich eflfect power spectrum (SZ). 
The two errors for the WL analysis by Refregier et al. (2002) are for statistical uncertainty and cosmic variance, respectively. 



cedure allowed us to have hints about the eflfect of treating 
the uncertainties in the L-T-M scaling relations in different 
ways. On the other hand, in order to treat errors adequately, 
one should be able to precisely assess their nature, and at 
the moment the split between scatter and overall normal- 
ization uncertainties is not well characterised, nor are these 
uncertainties precisely quantified. We notice, for example, 
that part of the previous disagreement between PSW and 
BOl was not in fact due to the adopted value of the M-T 
normalization T, (which would have given a discrepancy in 
the opposite direction), but rather to the statistical proce- 
dure used. 

For the XTF approach, the main scaling relation in- 
volved is the M-T one. Despite the fact that the agreement 
between observations and simulations has improved, mainly 
due to better modeling of the cluster physics in the simula- 
tions (e.g. Muanwong et al. 2002; Dave et al. 2002; Tornatore 
et al. 2002, in preparation), there is still a significant scatter 
in the quoted T, values. Should that be taken as the 'system- 
atic' error, it would certainly dominate the statistical error 
quoted by either observers or simulators (which is typically 
about 10 per cent). At the same time, observational mass de- 
terminations are probably affected by significant biases, due 
for instance to the necessary extrapolation from the (small) 
internal region observed to the larger one involved in the 
mass function. Other source include possible deviations of 
the surface brightness profile from simplified fitting mod- 
els and to the violation of the often made hydrostatic and 
isothermal assumptions. 

The same general comments apply to the XLF ap- 
proach. In this case, the fairly large number of clusters in- 
cluded in nearby and distant samples means statistical noise 



is unimportant for the resulting uncertainties. In fact the 
same complex physical processes which affect the ICM tem- 
perature are also very important when dealing with X-ray 
luminosity, which is highly sensitive to the local properties 
of the ICM. Therefore, more robust constraints on ag will re- 
quire substantial efforts to improve both the observational 
picture and the theoretical understanding of the relevant 
physical processes, which determine the X-ray properties of 
the intra-cluster medium. 

Because of the variety of ag values in the literature, 
we carried out careful tests of the impact of different £is- 
sumptions on the final result. On the theoretical side, it is 
important to rely on the most accurate mass function for 
the model being considered. In our analysis, we have used 
the .Jenkins ct al. (2001) mass function, and found that it 
gives closer results to pure Press-Schechter than to Sheth- 
Tormen for the mass range probed. Moreover, the use of the 
explicit as quantity instead of fitting formulae, as well as an 
accurate transfer function which properly includes baryonic 
suppression, is also necessary in order to reach the desired 
accuracy. While such effects are probably small with respect 
to the systematic uncertainties in the M—T—L relations, nev- 
ertheless they should be properly understood if the cluster 
abundance normalization is to be determined with high pre- 
cision. 

Bearing in mind all these caveats, a few firm conclusions 
can be drawn from our analysis. Firstly, current data on the 
XLF and XTF for nearby clusters give quite consistent re- 
sults on cosmological parameters, with erg — 0.8 for a flat 
Qm = 0.3 Universe. While one would be tempted to attach a 
~ 5 per cent statistical error (at la) to this determination, it 
is fair to say that at least twice as large an error, associated 
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with uncertainties in the M-T and, especiaUy, in the M- 
L relations, is more realistic. Secondly, under the same as- 
sumptions for the M-L conversion in the XLF analysis, dy- 
namical constraints on the {Q,m, cs) plane, derived from the 
evolution of the cluster population (e.g. BOl), are fully con- 
sistent with the geometrical constraints derived from nearby 
samples. 

The cluster abundance as a cosmological probe has been 
developing very fast in the past few years. New observa- 
tions have allowed an empirical determination of the scaling 
relations involved, while more refined simulations are help- 
ing to reconcile past discrepancies with observations. Fur- 
thermore, great effort has been devoted in understanding 
which modelling should be used in terms of transfer and 
mass function, and to the statistical procedures to compare 
observed and predicted number density of galaxy clusters. 
Furhcr improvements arc expected also on the observational 
side, thanks to the improved capability of X-ray satellites of 
the last generation to precisely measure cluster temperar 
tures. A cereful study over a statistically significant sam- 
ple of cluster will therefore allow a better calibration of the 
scaling relations and of their intrinsic scatter. At the same 
time, improvements in simulations, in terms of both reso- 
lution and physics included, will allow a reliable framework 
for the interpretation of the observational results. 

X-ray surveys, both in contiguous areas and based on 
cluster searches from Chandra and XMM archives, will cer- 
tainly help in improving the quantity and quality of cluster 
data. Future surveys, using missions like WFXT (Burrows 
et al. 1992) or DUET, as well as Sunyacv-Zcl'dovich sur- 
veys (sec Table 1 of Schulz & White 2002 for a list), will 
be extremely valuable for increasing the statistics of clus- 
ters at 2: > 1 by orders of magnitude, to enable evolutionary 
studies to be carried out. However, what is really needed is 
a robust set of mass estimates for a sample of clusters. An 
inter-comparison of multi-wavelength mass estimates, from 
X-ray data, velocity dispersions, gravitational lensing and 
Sunyaev-Zel'dovich measurements should enable systematic 
biases to be much better understood than they are at the 
moment. On the theoretical side, the challenge for the fu- 
ture will be to realize a new generation of cluster simula- 
tions, where the inclusion of the relevant physical processes 
should be coupled with a large enough dynamical range so 
as to accurately resolve the internal structure of the clusters 
while encompassing cosmological volumes. It will require a 
concerted effort from both the theoretical and observational 
communities, tackling the many remaining obstacles, in or- 
der to hone the cluster abundance into a precision cosmo- 
logical tool and to keep it competitive with rapid advances 
in the other methods. 
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